Measurements in the Levy quantum walk 
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We study the quantum walk subjected to measurements with a Levy waiting-time distribution. We 
find that the system has a sub-ballistic behavior instead of a diffusive one. We obtain an analytical 
expression for the exponent of the power law of the variance as a function of the characteristic 
parameter of the Levy distribution. 
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I. INTRODUCTION 

The development of the quantum walk (QW) in the 
context of quantum computation, as a generalization of 
the classical random walk, has attracted the attention of 
researchers from different fields. The fact that it is possi- 
ble to build and preserve quantum states experimentally 
has led the scientific community to think that quantum 
computers could be a reality in the near future. On the 
other hand, from a purely physical point of view, the 
study of quantum computation allows to analyze and ver- 
ify the principles of quantum theory. In this last frame 
the study of the QW subjected to different sources of 
decoherence is a topic that has been considered by sev- 
eral authors fl*]. In particular we have recently studied 
[^] the QW and the quantum kicked rotor in resonance 
subjected to noise with a Levy waiting-time distribution 
[s'l, finding that both systems have a sub-ballistic wave 
function spreading, as shown by the power-law tail of the 
standard deviation {a{t) ^ with 0.5 < c < 1), instead 
of the known ballistical growth {a{t) ^ t). This sub- 
ballistic behavior was also observed in the dynamics of 
both the quantum kicked rotor and the QW when 
these systems are subjected to an excitation that follows 
an aperiodic Fibonacci prescription. Other authors also 
investigated the kicked rotor subjected to noises with a 
Levy distribution Q and almost-periodic Fibonacci se- 
quence , showing that this decoherence never fully de- 
stroys the dynamical localization of the kicked rotor but 
leads to a sub-diffusion regime for a short time before 
localization appears. All these mentioned papers have in 
common that they work with quantum systems that have 
an anomalous behavior that was established numerically. 
There are no analytical results that explain in a general 
way why noises, with a power-law distribution or in a Fi- 
bonacci sequence, lead the system to a new non-diffusive 
behavior. Here we present a simple model that allows 
an analytical treatment to understand the sub-ballistic 
behavior. We hope that this may help to understand in 



a generical way how the frequency of the decoherence is 
the main factor in this unexpected dynamics. With this 
aim we investigate the QW when measurements are per- 
formed on the system with waiting times between them 
following a Levy power-law distribution. We show that 
this noise produces a change from ballistic to sub-ballistic 
behavior and we obtain analytically a relation between 
the exponent of the standard deviation and the charac- 
teristic parameter of Levy distribution. 

The paper is organized as follows. In the next section 
we develop the QW model with Levy noise, in the third 
section analytical results are obtained and in the last 
section we draw the conclusions. 



II. QUANTUM WALK AND MEASUREMENT 

The dynamics of the QW subjected to a series of mea- 
surements will be generated by a large sequence of two 
time-step unitary operators Uq and Ui as was done in a 
previous work |2|. But now C/q is the 'free' evolution of 
the QW and Ui is the operator that measures simultane- 
ously the position and the chirality of the QW. The time 
interval between two applications of the operator Ui is 
generated by a waiting-time distribution p{T), where T 
is a dimensionless integer time step. The detailed mech- 
anism to obtain the evolution is given in We take 
p{T) in accordance with the Levy distribution [1, [^, 
that includes a parameter a, with < a < 2. When 
a < 2 the second moment of p is infinite, when a = 2 
the Fourier transform of p is the Gaussian distribution 
and the second moment is finite. Then, this distribution 
has no characteristic size for the temporal jump, except 
in the Gaussian case. The absence of scale makes the 
Levy random walks scale-invariant fractals. This means 
that any classical trajectory has many scales but none in 
particular dominates the process. The most important 
characteristic of the Levy noise is the power-law shape of 
the tail, accordingly in this work we use the waiting-time 
distribution 
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To obtain the time interval T we sort a continuous vari- 
able t in agreement with eq. ffl and then we take the 
integer part Ti of this variable [2| ■ 

To obtain the operator C/q we develop in some detail 
the free QW model. The standard QW corresponds to 
a one-dimensional evolution of a quantum system (the 
walker) in a direction which depends on an additional 
degree of freedom, the chirality, with two possible states: 
'left' \L) or 'right' \R). Let us consider that the walker 
can move freely over a series of interconnected sites la- 
beled by an index n. In the classical random walk, a 
coin flip randomly selects the direction of the motion; in 
the QW the direction of the motion is selected by the 
chirality. At each time step a rotation (or, more gener- 
ally, a unitary transformation) of the chirality takes place 
and the walker moves according to its final chirality state. 
The global Hilbert space of the system is the tensor prod- 
uct Hs (Xi He where Hg is the Hilbert space associated to 
the motion on the line and He is the chirality Hilbert 
space. 

If one is only interested in the properties of the prob- 
ability distribution it suffices to consider unitary trans- 
formations which can be expressed in terms of a single 
real angular parameter 9 [ll|, [H, [l^. Let us call M_ 
(M+) the operators that move the walker one site to the 
left (right) on the line in Hs and let \L){L\ and \R){R\ be 
the chirality projector operators in He- Then we consider 
free evolution transformations of the form [ll[, 

Uo{e) = {Af_ ® \L){L\ + M+ ® \R) {R\} o {I ® K{e)} , 

(2) 

where K{d) = is an unitary operator acting on 

He, (Jy and cr^ being the standard Pauli matrices, and 
/ is the identity operator in H^. The unitary operator 
Uo{9) evolves the state \^{t)) by one time step, 

\^{t + i))^u,m'^{t)). (3) 

The wave- vector |^(i)) is expressed as the spinor 

(:;S)i..>. (4) 



that the two qubit states (l,i)/-\/2, and (1,— i)/\/2 are 
eigenstates of the measurement operator. This means 
that we project the chirality on the y direction using 
the ay Pauli operator. This form of measurement en- 
sures that the initial conditions after each measurement 
are equivalent to each other from the point of view of 
the probability distribution [14J. In this work we take 
6* = 7r/4 as in the usual Hadamard walk on the line. The 
probability distribution for the walker's position at time 
t is given by 



Pn{t) = |a„(t)|' + \bnit)\' 



(6) 



The effect of performing measurements on this system 
at time intervals Ti, with a Levy distribution, combined 
with unitary evolution is as follows: the standard de- 
viation has a ballistic growth during the free evolution, 
and when a measurement collapses the wave function the 
standard deviation starts again from zero. Then the stan- 
dard deviation has a zig-zag path. In Fig. [T] we present 
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where we have associated the upper (lower) component 
to the left (right) chirality, the states \n) are eigenstates 
of the position operator corresponding to the site n on 
the line. The unitary evolution for |^(t)), corresponding 
to eq. ^ can then be written as the map 



a„(t + 1) = a„+i(t) COS0 + bn+i{t) sinl 
bn{t + 1) ~ a„_i(t) sin0 — 6„_i(t) cos( 



(5) 



To build the measurement operator Ui we consider the 
case in which the walker starts from the position eigen- 
state |0), and with an initial qubit state (ao,6o) = 
(l,«)/-\/2. The operator Ui must describe the measure- 
ment of position and chirality simultaneously. The mea- 
surement of position is direct, but among the many ways 
to measure the chirality we choose to do it in such a way 



FIG. 1: The standard deviation for the QW as a function of 
dimensionless time in logarithmic scales. The parameters of 
the curves from top to bottom are: (1) a — 0.1 and c = 1.; 
(2) a = 0.5 and c = 0.99; (3) a = 1. and c = 0.92; (4) a = 1.5 
and c = 0.74 and (5) a = 2. and c = 0.51. 

the numerical calculation of the average standard devi- 
ation {a (t)) of the QW with measurements, calculated 
through a computer simulation of the time evolution of 
an ensemble of 2 x 10^ stochastic trajectories for each 
value of the parameter a. This figure shows for different 
values of a that the behavior is diffusive (^ t) for times 
t < 10 and follows a power law ~ for times t > 10. One 
would expect that the large degree of decoherence intro- 
duced by measurements led to a diffusive behavior for all 
times [2!|; instead the system changes to an unexpected 
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sub-ballistic regime. Then a determines the degree of 
diffusivity of the system for large waiting times, that is 
ballistic for a = 0, sub-ballistic for < a < 2 and dif- 
fusive for a — 2. Fig. [5] shows, in full line, the exponent 
c of the power law as a function of the Levy parameter 
a. Again the calculation has been made with an ensem- 
ble of 2 X 10^ stochastic trajectories for each value of a. 
The passage from the ballistic behavior for a = to the 
diffusive behavior for a = 2 is clearly observed. In this 
figure the result of the analytical calculation of the next 
section is also presented in dashed line, the comparison 
of the curves proving the coherence of both treatments. 

III. THEORETICAL MODEL 

Different mechanisms of unitary noise may drive a 
system from a quantum behavior at short times to a 
classical-like one, at longer times. It is clear that the 
quadratic growth in time of the variance of the QW is 
a direct consequence of the coherence of the quantum 
evolution [15]. In this section, we develop an analytical 
treatment to understand why the Levy decoherence does 
not fully break the ballistic behavior. 

In our previous work [l^ we investigated the QW on 
the line when decoherence was introduced through simul- 
taneous measurements of the chirality and position. In 
that work it was proved that the QW shows a diffusive 
behavior when measurements are made at periodic times 
or with a Gaussian distribution. Now, we use the Levy 
distribution but some of the result obtained in [l6l | can 
be used. For the sake of clarity we reproduce briefly the 
main steps to obtain the dynamical equation of the vari- 
ance. Let us suppose that the wave-function is measured 
at the time t, then it evolves according to the unitary 
map ([5]) during a time interval T, and again at this last 
time t + T, a new measurement is performed. The prob- 
ability that the wave-function collapses in the eigenstate 
|n) due to the position measured, after a time T, is 

qn = P,,{T). (7) 

These spatial distributions (/„ depend on the initial qubit 
state and the time interval T, they will play the role 
of transition probabilities for the global evolution. The 
mechanism used to perform measurements of position 
and chirality assures that these distributions will repeat 
themselves around the new position, because the ini- 
tial chirality ((l,i)/\/2or (1,— i)/\/2) produces the same 
spatial distributions g„ and their value only depends on 
the size of T. Then it is straightforward to build the 
probability distribution P„ at the new time i -I- T as a 
convolution between this distribution at the time t with 
the conditional probability qn(T); this takes the form of 
the following 'master equation' 

n+T 

Pn{t + T)= ^ g„-,P,(t), (8) 

3=n~T 



where Qn-j are the transition probabilities from site j 
to site n, defined in eq. ([7|) and the sum is extended 
between j = n — T and j = n + T, because T is also 
used as the number of applications of the quantum map 
eq. ([5]). Remember that for each time step the walker 
moves one spatial step to the right and left. Using eq. ([5]) 
we calculate the first moment mi{t) = ^jPj{t) and the 
second moment m2{t) = J^J^Pji^) to obtain 

mi{t + T) =mi{t) +mig{T) (9) 
m2{t + T)= m2{t) + 2mi{t)miq(T) + m2q(T) (10) 

n=T n=T 

where mig(r) = ^ ng„ and m2q{T) — ^ ri^qn are 

n=-T n=~T 

the first and second moments of the unitary evolution 
between measurements. Therefore the global variance 
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FIG. 2: the exponent c of the power law of the standard 
deviation for the QW as a function of the parameter a. The 
full line correspond to the numerical result and the dashed 
line to the analytical one 

cr^(t) = 7712 (t) — (i) verifies the following equation for 
the process, obtained for first time in |16|] 

a^{t + T) = a^{t) + al{T), (11) 

where <Tg{T) = m2q{T) — m\q{T) is the variance asso- 
ciated to the unitary evolution between measurements. 
These results can be used for random time intervals T 
between consecutive measurements. The variance ct^ de- 
pends very weakly on the qubit 's initial conditions and 
it increases quadratically with time [Tlj 

cjl{T) = kT\ (12) 

where T ^ 1 and fc is a constant determined by the 
initial conditions. We now calculate the average of eq. 
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(fTTjl taking integer time intervals Ti between consecutive 
measurements sorted according to the Levy distribution 
eq. (dl) 



{a\t + T,)) = {a^t)) + k{T^), 



(13) 



where {f(t)) = J f {x)p{x)dx. From the previous nu- 



merical calculation, we know that the variance grows as 
a power law, then (cr-^ (t)) cx t'^'^ and (^a'^{t + Ti)) cx 

((i + r,)2'=) « + 2c^) for large t. Substituting 

these expressions in eq. the following result for the 
exponent c is obtained 



(14) 



valid for a large t. The first and the second moments of 
the waiting time of our Levy distribution are 
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(15) 



J-2-Q 



2-a 



(16) 



Therefore, in the case when t — > oo the exponent c is 



1 , if ^ a ^ 1 
^~ i(3-a), if 1 a 2. 



(17) 



This result is in accordance with the numerical result ob- 
tained in the previous section, as can be seen in Fig. [2l It 
is important to remark that eq. (jl7p gives the analytical 
dependence of the exponent of the power law for (ct^) on 
the parameter a. 



IV. CONCLUSION 

Several systems have been proposed as candidates to 
implement the QW model, they include atoms trapped in 
optical lattices (itI [Tsj , cavity quantum electrodynamics 
IS I and nuclear magnetic resonance in solid substrates 
2C , im . All these proposed implementations face the ob- 
stacle of decoherence due to environmental noise and im- 
perfections. Thus the study of the QW subjected to dif- 
ferent types of noise may be important in future techni- 
cal applications. Here we showed that the QW subjected 
to measurements with a Levy waiting-time distribution 
does not break completely the coherence in the dynam- 
ics, but produces a sub-ballistic behavior in the system, 
as an intermediate situation between the ballistic and the 
diffusive behavior. Note that as Gaussian noise is a par- 
ticular case of the Levy noise, our study is open to wider 
experimental situations. We studied this behavior nu- 
merically and we obtained also an analytical expression 
for the exponent of the power law of the variance as a 
function of the characteristic parameter of the Levy dis- 
tribution. Measurement is one of the strongest possible 
decoherences on any quantum system, but even in this 
case the coherence of the QW treated here is not fully 
lost . This fact shows that the most important ingredi- 
ent for the sub-ballistic behavior is the temporal sequence 
of the perturbation, not its size or type. This extreme 
and simple model that can be treated analytically allows 
us to understand why similar models 0, 0, IE] (that were 
studied numerically) show also a power-law behavior. 
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